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Abstract 


Continuous  phase  modulation  (CPM)  is  a  well-established  signaling  technology  that 
attains  spectral  efficiency  by  smoothing  the  phase  signal.  The  spectral  efficiency  of 
CPM  makes  it  a  natural  candidate  for  communications  where  bandwidth  is  limited, 
such  as  the  non-SATCOM  channels  in  the  littoral  environments.  However,  CPM 
receivers  require  synchronization  with  respect  to  data,  timing,  phase,  and  frequency. 
Research  on  various  CPM  synchronizers  has  been  undertaken  at  Space  and  Naval 
Warfare  Systems  Command  (SPAWAR).  This  report  sets  out  a  simulation  framework 
that  allows  these  synchronizers  to  be  swapped  in  and  out  of  a  CPM  receiver  model 
to  assess  performance  in  a  credible  end-to-end  RF  simulation. 
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1 

The  CPM  Simulation  Framework 


Continuous  phase  modulation  (CPM)  is  a  mature  signaling  technology  that  attains 
spectral  efficiency  by  smoothing  the  phase  signal.  The  tradeoff  is  that  such  a  “slip¬ 
pery”  signal  can  be  hard  to  decode.  From  Sklar  [34,  page  631]: 

The  bandwidth  efficency  of  CPM  is  obtained  by  increasing  the  smoothness 
of  the  waveforms  in  the  time  domain.  . . .  However,  this  smoothness  in  the 
time  domain  also  tends  to  eliminate  the  symbol  transition  features  upon 
which  many  symbol  synchronization  schemes  depend. 


In  particular,  CPM  signal  processing  requires  that  the  signal  be  synchronized  with 
respect  to 


•  timing 

•  frequency 

•  phase 


before  decoding  is  attempted.  The  (complex  baseband)  CPM  signal  has  the  form  [34, 
Section  10.2.3]: 

sB(t)  =  exp  (JMt)), 


where  the  phase 

OO 

=  2tt  53  hakq(t.  -  kT ) 

k= — oo 

carries  the  information  in  the  ak  s.  These  symbols  modulate  the  continuous  phase 
pulse  q(t).  Figure  1.1  shows  a  specific  phase  pulse  to  illustrate  the  general  shape  of 
the  phase  pulse.  Because  the  q(t)  is  continuous,  the  phase  VleW  is  continuous.  The 
smoothness  of  ip  sit)  and  the  modulation  index  h  determine  the  bandwidth  of  the 
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Figure  1.1:  Phase  pulse  1REC;  T  =  1. 

CPM  signal.  The  received  CPM  signal  is  corrupted  by  additive  noise  and  demodula¬ 
tion  errors: 

rB(t)  =  ejA^ej2nAfctsB{t  -  r)  +  gB(t). 

The  synchronizers  must  compensate  for  the  phase  A<f>,  frequency  A /c,  and  timing  r 
offsets. 

Multi-/?,  CPM  originated  in  1978  paper  of  Anderson  &  Taylor  [5].  Multi-/?,  CPM 
generalizes  single-/?  as 

OO 

^B(t)  =  27 t  hk(*kq{t  -  kT ). 

k=—OG 

Consequently,  multi-/?,  also  requires  synchronization  on  the  modulation  indices.  Com¬ 
pared  to  single-/?,  multi-/?  codes  can  reduce  the  bandwidth  [21]  and  increase  the 
probability  of  detection  [7,  Section  3.2.3]. 

Synchronizers  are  required  to  process  the  CPM  signals.  The  idea  is  that  phase  Ac/?, 
frequency  A/,  and  timing  r  are  compensated  so  that  a  good  decoding  of  the  phase  is 
possible.  Typically,  each  synchronizer  is  implemented  with  a  Vit.erbi  decoder.  Several 
decoders  may  be  used  leading  to  complex  CPM  receivers  [31].  Figure  1.2  shows  such 
a  layout  for  a  multi-/?,  CPM  receiver  [36],  [21].  that  handles  phase,  symbol  timing, 
and  superbaud  timing  in  addition  to  the  symbol  decoding. 

As  expected,  the  literature  on  synchronization  is  large  with  papers  filling  in  the 
matrix  whose  elements  are  indexed  by  synchronizers,  waveforms,  multi-/?,  codes,  re¬ 
ceivers,  channels,  and  applications.  However,  we  are  not  interested  in  the  particulars 
of  a  synchronizer  but  care  only  that  it  improve  the  end-to-end  performance  of  the 
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Figure  1.2:  Classical  multi-/?.  CPM  receiver. 

CPM  receiver.  For  example,  suppose  a  timing  synchronizer  could  be  tweaked  to  get 
a  10%  reduction  in  timing  error.  Is  this  significant?  Only  if  this  reduction  in  timing 
error  improves  performance  of  the  CPM  radio.  So  the  goal  of  this  report  is  to  develop 
a  simple  CPM  radio  that  allows  us  to  swap  out  the  various  synchronizers. 

To  meet  this  goal,  we  developed  of  series  of  simple  CPM  receivers  CMPO,  CPM1, 
and  CPM2  that  permit  us  to  “plug-n’-play”  the  synchronizers.  Section  2  starts  by 
reviewing  the  basics  of  the  multi-/?  signals  and  sets  out  the  very  simple  CPMO  re¬ 
ceiver.  Section  3  reviews  selected  synchronizers  and  reports  on  their  performance  at 
estimating  the  demodulation  errors.  Section  4  generalizes  CPMO  to  handle  timing 
and  phase  in  CPM1.  Timing,  phase,  and  frequency  are  handled  in  CMP2.  Section  5 
assess  the  performance  of  multi-/?  CPM  in  multipath.  The  channel  considered  is  one 
of  the  military  UHF  models  for  a  littoral  zone.  For  the  particular  implementation  of 
the  multi-/?,  the  multipath  is  essentially  a  noise  source  so  communication  is  degraded. 
There  exist  several  solutions  (i.e.,  demodulate  both  the  signal  and  channel).  A  review 
of  these  solutions  and  associated  research  topics  concludes  this  report. 
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2 

Multi-/?  Signals 


This  section  sets  the  notation  and  nomenclature  for  processing  the  multi-/?,  continuous 
phase  modulation  (CPM)  signals.  Sections  2.1  and  2.2  review  the  transmitted  and 
received  multi-/?,  signal  at  baseband.  The  received  signal  is  distorted  by  demodulation 
errors  and  corrupted  by  additive  noise.  The  additive  noise  is  set  by  the  signal-to-noise 
ratio  (SNR).  However,  receiver  performance  is  typically  measured  by  plotting  the  error 
rate  against  Eb/N0  rather  than  the  SNR.  To  map  the  SNR  to  Eb/N0  requires  that 
we  know  the  spectra  of  the  multi-/?  signals.  The  multi-/?  codes  determine  the  multi-/? 
spectra  as  illustrated  in  Section  2.3.  These  sections  complete  the  simulation  of  the 
received  multi-/?  signal. 

Processing  the  multi-/?  signal  requires  synchronizing  in  spite  of  the  demodulation 
errors.  The  processing  requires  compensating  for  the  demodulation  errors  and  can 
lead  to  complex  CPM  receivers.  Because  we  want  the  end-to-end  performance  of  the 
synchronizers,  Section  2.4  offers  a  simple  CPM  receiver.  The  simplicity  of  this  receiver 
allows  us  to  swap  out  the  various  synchronizers  to  get  the  end-to-end  comparisons. 


2.1  The  Transmitted  Multi-/i  Signal 

The  baseband  multi-/?,  signal  sB(t)  has  the  form  [10]: 

OO 

sB(t)  =  exp(jV>js(Z));  =  27 t  hyk]akq(t  -  k.T).  (2.1) 

k=— oo 

The  modulation  indices  {/?0,  /?i, . . . ,  /?//  i }  typically  cycle  as 

[A;]  =  mod(A:,  H). 

The  symbol  sequence  {ct*.}  consists  of  independent  elements  drawn  from  the  alphabet 

ak  £  {±1,  ±3, . . . ,  ±(M  -  1)}, 
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where  M  is  even.  Each  symbol  is  transmitted  over  the  symbol  period  T.  The  pnlse 
q(t)  is  known  as  the  phase  pulse  [24],  or  the  pulse-shaping  pulse  [10],  or  the  modulation 
phase  response  [34],  The  phase  pulse  q(t. )  is  the  integral  of  the  frequency  pulse  q'(t) 
[24]: 

?(*)  =  /  (l'{r)dr. 

The  frequency  pulse  is  supported  on  the  interval  [0,  LT]  and  normalized  so  that 


Q(t)  = 


0  t  <  0 

1/2  t  >  LT  ’ 


where  L  is  called  the  correlation  length  [34,  page  633].  When  L  —  1,  the  signaling  is 
called  full  response.  [34,  page  633].  If  L  >  1,  the  phase  pulses  overlap  and  the  signaling 
is  called  partial  response. 


2.2  The  Received  Multi-/i  Signal 

The  received  baseband  signal  rB(t)  is  corrupted  by  additive  noise  and  demodulation 
errors: 

rB(t )  =  e,**ej2"AfetsB(t  -  r)  +  gB(t). 

•  A  4>  :=  phase  error 

•  A/e  :—  frequency  error 

•  t  :=  timing  error 

Jitters  in  the  phase,  the  frequency,  or  the  timing  are  not  modeled.  The  additive  noise 
{9b(1)}  is  zero-mean,  complex- valued,  wide-sense  stationary  bandlimited  Gaussian 
noise  with  variance  determined  by  the  signal-to-noise  ratio: 

C'NT'n  \  _ \ _ 

Var  [gB]  E[\gBm' 

We  assume  the  received  signal  is  sampled  at  Nsps  samples  per  symbol.  That  is,  the 
symbol  period  T  is  a  multiple  of  the  sample  period  Ts : 

T  =  NspsTs. 

The  digital  received  signal  is  the  vector 

x  =  [rB(kTs)}  k  =  0,1, . . .  Nt  —  1. 

In  this  setup,  the  sample  rate  fs  is  a  multiple  of  the  CPM  “bandwidth” 

N 

n  1  y  sps 
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so  that  aliasing  of  sB(t)  may  be  ignored.  If  the  received  signal  is  observed  over  Nsym 
symbols,  a  total  of  N,  samples  is  obtained.  The  total  observation  period  is 

N,TS  =  NsymNspsTs  =  NsymT. 

Table  2.1  summarizes  the  sampling  notation. 

Table  2.1:  Sampling  notation. 


NSym  Number  of  symbols  observed 
Nt  Number  of  time  samples 
iVSpS  Number  of  samples  per  symbol 
T  Symbol  period 
Ts  Sampling  period 


The  simulations  take  the  noise  {gs{t)}  bandlimited  to  the  frequency  interval 
[— /s/2,  fs/2].  Consequently,  the  sampled  noise  {g£>(kTs)}  is  a  time  series  that  is 
zero  mean,  complex- valued,  IID  Gaussian  with  variance  Var  [(/#].  As  the  sample  rate 
fs  increases  and  the  SNR  is  held  constant,  the  noise  in  the  CPM  band  decreases. 

The  ubiquitous  E\,/Nq  is  a  normalized  figure  of  merit  for  digital  radios  [34,  Section 
3.1.5].  The  SNR  is  mapped  to  E^/No  as  [34,  Eq.  3.30]: 

f  =  SNRx  (2.2) 

iVo  Kb 

where  W  and  Ri,  denote  the  bandwidth  (Hz)  and  bit  rate  (bits/sec),  respectively. 
Computing  the  bit  rate  R\,  is  straightforward.  With  the  alphabet  {±1,  ±3,...  ± 
(M  —  1)},  there  are  log2(M)  bits  per  symbol  so  the  bit  rate  is  [34,  Eq.  9.14]: 

%= (Hz). 


Computing  the  bandwidth  is  not  as  straightforward  because  “bandwidth  does  not 
have  a  unique  definition”  [37],  [34,  Figure  1.20].  Sundberg  uses  the  99%  bandwidth 
/g9  [37],  as  do  Ho  &  Mclane  [14]: 


0.99 


ht  r-Wf 


Here  Pas(f )  denotes  the  power  spectrum  of  the  CPM  signal.  We  adopt  /g g  as  the 
bandwidth  of  the  multi-/?  CPM  signal.  Consequently,  multi-h  spectrum  must  be 
computed  for  each  modulation  index. 
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2.3  Multi-/i  Spectra 

How  the  modulation  indices  determine  the  multi-/?,  spectra  is  a  fascinating  topic  that 
has  led  to  a  variety  of  computational  methods  [14],  [7,  Chapter  4],  [37],  The  direct 
method  simply  estimates  the  multi-/?,  spectrum  from  a  simulation  [44] .  The  following 
figures  plot  these  spectral  estimates  for  selected  multi-/?  indices.  The  simulations  use 
symbol  and  sampling  periods  set  to  the  nominal  values: 

T=  1,  Ts  =  10. 


The  number  of  symbols  is 

fVsym  =  1000. 

Estimates  of  the  95%  confidence  intervals  bound  the  spectrum.  (These  confidence 
intervals  are  computed  pointwise  under  an  asymptotic  normality  assumption  [20].) 
Also  reported  on  each  plot  is  an  estimate  of  the  99%  bandwidth.  This  bandwidth  is 
used  to  map  the  SNR  to  Eb/No- 
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CPM  Spectra  for  /?=[8/16  14/16]  @  95%  confidence  intervals 
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Figure  2.2:  Spectral  estimates  of  [/?o  hi]  =  [8  14]/16. 


CPM  Spectra  for  /;=[  12/16  15/16]  @  95%  confidence  intervals 


frequency  (Hz):  N= 10000;  A/sym=1000 


Figure  2.3:  Spectral  estimates  of  [h0  /ij]  =  [12  1 5] / 1 6 . 


10 


2.4  A  Multi-/i  Receiver 


Typical  CPM  radios  are  baroque  designs  [21],  [31]  with  one  or  more  Viterbi  decoders 
devoted  to  each  of  the  following  tasks:  symbol  timing  [6],  superbaud  timing  [10], 
frequency  estimation  [4],  phase  estimation  [23] — not  to  mention  the  Viterbi  decoder 
for  the  data  [19].  Instead  of  following  these  well-worn  designs,  we  offer  a  simple  CPM 
receiver  whose  transparent  design  lets  us  swap  the  synchronizers  in  and  out,  shows  the 
effect  of  multipath,  adapts  to  the  Walsh  receivers,  and  reveals  particular  difficulties 
of  CPM  processing. 

We  develop  a  series  of  CPM  receivers,  CMPO,  CPM1,  CPM2,  . .  .of  increasing 
capability.  The  starting  point  is  a  very  simple  CPM  receiver.  CPMO  assumes  perfect 
synchronization  in  additive  noise: 

•  No  timing  errors 

•  No  phase  errors 

•  No  frequency  errors 

CPMO  assumes  that  Nsym  symbols  have  been  transmitted 

Afsym 

sB(t.)  =  ));  ipB(t)  =  2tt  Y  h{k]akq(t  -  kT ) 

k= 1 

and  that  the  received  baseband  signal  is  corrupted  only  by  additive  noise: 

rB{t)  =  sB{t)+gB(t). 

The  problem  is  to  recover,  and  decode  ipB(t)  from  rB(t).  CPMO  simply  unwraps  the 
phase,  estimates  the  instantaneous  frequency,  and  decodes  by  least  squares: 

CPMO-1  Estimate  the  phase  of  sB(t)  from  rB{t). 

CPMO-2  Estimate  the  instantaneous  frequency  from  the  estimated  phase. 

CPMO-3  Estimate  the  alphabet  from  the  instantaneous  frequency. 

Introduce  the  piecewise-linear  function 

-Nsym 

m  ■■=  2tt  £  fkq(t  -  kT)2T. 

k= 1 

With  L  =  1,  the  first  derivative  is 

OO 

£'(t)  :=  2tt  Y,  fklv>,T\(t  -  kT). 

k~ — oc 

Let.  jC  denote  collection  of  all  these  piecewise  linear  functions — C  is  then  a  linear  space 
of  dimension  Nsym.  The  CPMO  design  is  detailed  as  follows: 
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CPMO-1  Estimate  the  phase  of  sB(t)  by  phase  unwrapping  the  received  signal  [20]: 

:=  unwrap(angl e(rB(t))). 

CPMO-2  Find  the  best  piecewise  linear  approximation  to  the  estimated  phase: 

1  ■=  argmin{||^s  -  ^||2  :  i  €  £}. 

CPMO-3  Decode  the  phase.  With 


■/Vsym 

((f)  :=  2?r  y  Mt  -  kT)2T 

k= 1 

decode  by  the  best  fit,  in  each  symbol  interval: 

ak  :=  argmin{|/fc//?.[fc]  -  a\  :  a  =  ±1,  ±3} 

The  virtue  of  this  CPM  receiver  is  its  relative  simplicity — we  can  see  exactly  where 
the  synchronizers  fit.  This  approach  is  similar  to  the  polynomial-phase  signaling 
developed  in  the  mid-1990s  [8],  [30],  [15]. 

Figure  2.4  displays  the  performance  of  CPMO  for  selected  indices.  These  simula¬ 
tions  use  the  symbol  and  sampling  periods  set  to  the  nominal  values: 

T  =  1,  Ts  =  10. 

The  figure  plots  the  probability  of  symbol  error  Pe  as  a  function  of  SNR.  The  figure 
shows  a  considerable  difference  in  performance  as  a  function  of  the  indices.  However, 
most,  of  the  performance  plots  in  the  literature  compare  probability  of  bit  error  PB 
to  Eb/N0.  Thus,  we  need  to  convert  this  figure  to  make  comparisons  with  the  CPM 
literature. 

Figure  2.5  displays  the  performance  of  CPMO  as  a  function  of  Eb/Ng .  Equation  2.2 
makes  the  conversion  from  the  SNR  of  Figure  2.4  to  Eb/Ng  but  requires  the  bandwidth 
W  of  the  multi-/?  signal.  In  Section  2.2,  bandwidth  was  defined  as  W  =  fgg.  The 
spectral  plots  of  Section  2.3  show  how  the  bandwidth  is  controlled  by  the  multi-h 
indices  and  reports  on  our  estimate  of  the  fgg  bandwidth.  It  is  this  estimate  of  fgg 
that  maps  the  SNR  to  Eh/N0.  For  comparison,  two  multi- h  CPM  receivers,  also 
operating  with  perfect  synchronization,  are  also  plotted  on  Figure  2.5.  These  CPM 
receivers  from  the  literature  are  measured  using  the  bit.  error  rate  PB.  Assuming  the 
bit  errors  are  independent,  the  symbol  error  links  to  the  bit  error  rate  as: 

Pe  =  1  —  Prob(  no  bit  errors  ) 

=  1  -  Prob(l  -  PB)Nb'l% 


12 


10 


-1 

10 


-2 

10 


cC  10' 


10 

io'5 


10 

01234  5  6789  10 

SNR  (dB):  N=  100 

'  sym 

Figure  2.4:  SNR  performance  of  CPMO — perfect  synchronization. 

where  iVbits  is  the  number  of  bits  per  symbol.  This  is  the  formula  we  use  to  map 
between  Pb  and  Pe.  Figure  2.5  shows  that  the  performance  of  Premji  &  Taylor’s 
CPM  receiver  [31]  is  comparable  to  our  simple  CPMO  receiver.  The  extra  complexity 
of  Mazur  &  Taylor’s  receiver  [21]  pays  a  dividend  of  better  performance. 

We  close  this  section  with  one  way  to  rank  the  multi-h  indices.  This  ranking 
follows  from  the  Walsh  receivers  [38],  [39].  Assume  full  response  (L  —  1).  Each 
symbol  and  index  determines  a  signal  over  a  symbol  interval: 

s( h;  k,  m\  t)  =  exp(j2namh[k]q(t));  ( t  e  [0,  T]). 

For  matched  filtering,  ideal  performance  occurs  if  the  collection  of  signals  {.s(h;  k,  m)} 
is  orthogonal: 

r t _ 

0  =  /  s(h;  k,  m\  t)s(h;  k' ,  m'\  t)dt  (k,m)  ^  (k',m'). 

Jo 

For  sampled  CPM  signals,  let  s(h;  k,  m ;  t)  denote  vector  of  the  samples  of  the  ( k ,  m) 
signal: 

s(h;  k,  m;  t)  =  [s(h;  k ,  m;  nTs)  :  n  =  0, 1, . . . ,  Nsps  -  1]. 


CMP0:  Probability  of  Symbol  Error 
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CMPO:  Probability  of  Symbol  Error 


o 


Figure  2.5:  Eb/N0  performance  of  CPMO — perfect  synchronization;  o=Premji  &  Tay¬ 
lor  [31],  o=Mazur  &  Taylor  [21]. 

A  measure  of  the  orthogonality  is  determined  by  the  singidar-value  decompositions 
(SVD)  of  the  signal  matrix 

[s(h;  :  k  =  0, 1  ;m  =  1,2, 3,4]. 

Here  the  signal  matrix  is  made  specific  to  our  setup.  With  M  =  4  and  using  2-h,  there 
are  8  signal  vectors.  With  Nsps  =  10,  the  signal  matrix  is  10  x  8  and  has  singular  values 
s(l)  >  s(2)  >  . . .  >  s(8)  >  0.  If  the  8  signal  vectors  were  orthogonal,  s(8)/s(l)  =  1. 
If  the  signal  matrix  loses  orthogonality,  this  ratio  decreases.  If  the  ratio  decreases  to 
zero,  the  signal  matrix  is  rank  deficient.  Figure  2.6  sorts  the  s(8)/s(l)  ratio  of  the 
signal  matrices  as  a  function  of  the  multi-/;  indices.  The  most  orthogonal  index  is 
[12  15]/16  that  corresponds  to  the  maximal  point  on  the  curve.  Table  2.2  reports  the 
11  most  orthogonal  indices. 

However,  ranking  by  orthogonality  does  not  account  for  bandwidth.  The  spectral 
plots  indicate  that  more  orthogonality  costs  more  bandwidth.  Figure  2.7  makes  this 
observation  concrete  by  plotting  bandwidth  and  the  orthogonality  for  each  [h0  hi]. 
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Orthgonality  of  the  CPM  signals 


index:  optimal  [h0,hj]=[  12  15j/16 

Figure  2.6:  Measuring  the  orthogonality  of  the  [??i  712]/ 16  indices. 

A  few  points  are  listed  with  their  indices.  For  example,  this  plot  shows  that  [ho  hi]  = 
[4  5]/16  is  near  the  “boundary”.  As  such,  it  is  one  of  the  many  “optimal”  or  Pareto 
points.  This  plot  also  suggest  other  candidate  indices.  For  example,  [5  7]/16  achieves 
almost  10  times  the  orthogonality  over  [4  5]/16  with  only  a  modest  increase  in  band¬ 
width.  Thus,  we  expect  the  larger  bandwidth  lets  [5  7]/16  outperform  [4  5]/16  on  the 
SNR  plot.  Moreover,  we  suspect  that  this  better  performance  is  not  normalized  away 
by  a  much  larger  bandwidth  so  expect  [5  7]/ 16  outperforms  [4  5]/16  on  the  E^/No 
plot. 

The  following  figures  verify  these  expectations.  Figure  2.8  adds  [5  7]/16  to  the 
SNR  performance  plot  and  shows  it  lies  midway  between  the  SNR-optimal  curves 

Table  2.2:  Optimal  [h0  hi]  indices  ranked  by  orthogonality  (left-to-right,  corresponds 
to  least-to-most  orthogonality). 

h0  x  16  8  9  8  12  10  11  9  10  13  11  12 

hi  x  16  14  14  14  14  14  14  15  15  15  15  15 
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Muhi-/j  CPM  tradeoff:  Bandwidth  versus  Orthogonality 


0.6  0.8  1  1.2  1.4  1.6  1.8  2  2.2 

99%  bandwidth  (Hz) 


Figure  2.7:  Bandwidth  and  Orthogonality  tradeoffs  for  [ri\  n2]/16. 

and  the  [4  5]/16  curve.  Figure  2.9  adds  [5  7]/16  to  the  Eb/N0  performance  plot. 
Compared  to  [4  5]/16,  the  slight  increase  in  bandwidth  is  offset  by  the  improved 
detection.  However,  all  these  error  plots  assume  perfect  synchronization^the  real 
measure  of  any  CPM  receiver  is  how  well  it  can  extract  the  phase  information  when  the 
received  signal  is  not  synchronized.  Accordingly,  Section  3  is  devoted  to  synchronizers. 
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3 

Synchronizers 


The  received  baseband  CPM  signal  rB(t)  is  modeled  as  the  transmitted  baseband 
CPM  signal  sB(t)  distorted  by  the  demodulation  errors — phase  A(f>,  frequency  A/c, 
and  timing  r — and  corrupted  by  additive  noise: 

rB(t)  =  -  r)  +  gB(t). 

A  CPM  synchronizer  estimates  (A </>,  A fc,  r)  to  compensate  for  these  demodulation 
errors.  A  multi-/?  CPM  synchronizer  not  only  estimates  the  delay — the  superbaud 
timing — but  also  the  location  of  the  multi-/?  indices — the  symbol  timing.  The  most 
general  synchronizer  jointly  estimates  (A0,  A/c,  r).  Because  of  mathematical  diffi¬ 
culties,  the  parameters  are  also  estimated  separately:  phase  and  data  [19];  phase  and 
timing  [23];  frequency  and  data  [24], 

Synchronizers  have  been  developed  either  ad  hoc  or  a  formal  approach  based  on 
estimation  theory.  Both  approaches  have  proven  useful  and,  in  many  cases,  the  ad 
hoc  synchronizers  fit  within  an  estimation  theory  framework.  The  formal  approach  is 
often  based  upon  the  maximum,  likelihood  criterion  of  optimality — maximize  a  prob¬ 
ability  density  function  for  the  received  samples  conditioned  upon  the  parameters  to 
be  estimated  [41],  Because  of  compromises  in  both  the  analysis  and  receiver  imple¬ 
mentation,  an  actual  synchronizer  may  only  be  an  approximation  to  the  maximum 
likelihood  solution  and  consequently  suboptimal. 

Synchronizers  are  statistical  estimators.  A  measure  of  “goodness”  of  an  estimator 
is  its  variance,  usually  compared  to  the  Cramer-Rao  lower  bounds  [22].  We  care  only 
how  well  a  synchronizer  “works”  as  part  of  a  CPM  receiver  to  improve  end-to-end 
performance.  Consequently,  the  synchronizers  are  developed  to  be  swapped  in  and 
out  of  the  simple  CPM  receiver  of  Section  2. 

The  following  sections  develop  a  variety  of  synchronizers.  Section  3.1  reviews 
an  Italian  approach  to  superbaud  and  symbol  timing.  Section  3.2  covers  a  timing 
synchronizer  specialized  t.o  2-/?,.  Section  3.3  reviews  a  frequency  synchronizer.  The 
final  section  reviews  the  current  literature  on  synchronizers  and  points  out  more 
advanced  synchronizers  for  future  extensions. 
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3.1  Timing  Italian  Style 

The  Italian  school  claims  [10] :  “coherent  detection  of  multi-/?  CPM  requires  knowledge 
of  the  beginning  of  each  index  cycle.”  Their  contribution  is  the  representation  of  the 
delay  as 

r  =  (i 7  +  e)T. 

•  Superbaud  timing  parameter :  77  =  0, 1, ...  ,H  —  1. 

•  Symbol  timing  parameter.  —1/2  <  e  <  1/2. 

The  delay  r  is  estimated  by  estimating  the  symbol  timing  parameter  e  and  then  using 
e  t.o  estimate  rj. 

3.1.1  Symbol  Timing  Estimation 

An  estimator  e  for  symbol  timing  is  pieced  together  from  [10]  and  [6].  The  maximum 
likelihood  function  is  the  quadratic  form  [10]: 


A4(x|e)  =  xH[F(k]  -  k2)Ts ,  k,Ts  -  eT)]x, 


where  F  is  given  by  [1,  Eq.  3.2]: 

Ft  At  t)  =  TT  1  sm(2irMh[k]g(At,t  -  kT)) 

'  Jio  M  sm(2irh[k]q(M,  t  -  kT ))  ’ 

q(At.,t)  =  q(t)  -  q(t  -  At). 

The  simulations  show  that  A4(x|e)  tends  to  exhibit  a  single  maximum  so  that 

e  =  argmax{A4(x|e)  :  e  6  [—1/2, 1/2)} 

can  be  referred  to  as  “the”  maximum  likelihood  estimate. 

Figure  3.1  illustrates  the  performance  of  A4(x|e)  as  the  number  of  symbols  in  x  is 
increased.  For  comparison,  each  estimate  is  scaled  to  the  unit  interval: 

A4(x|e)  —  min{A4(x[e)} 

4  max{A4(x|e)}  —  min{A4(x|e)} ’ 

The  low-noise  simulation  shows  5  to  10  symbols  barely  work.  To  register  a  clean 
maximum,  we  see  that  20  to  30  symbols  are  needed.  This  “need”  brings  us  to  an 
excellent  research  topic: 

How  m.any  symbols  are  needed? 
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Maximum  Likelihood  Estimate  for  Symbol  l  iming  T-0 


e:  /H 4/16  5/16];  W  =10;  SNR-20  dB 


Figure  3.1:  A4(x|e)  at  r  =  0  using  x  ~  5,  10,  20,  and  30  symbols. 

This  question  offers  many  avenues  of  exploration.  We  can  ask  how  the  timing  error 
behaves  as  a  function  of  the  number  of  symbols.  We  also  ask  what  precision  in  e  is 
needed  to  synchronize.  Thus,  we  are  focused  less  on  c  and  more  on  its  performance 
in  a  CPM  receiver. 

To  show  the  effect  of  increasing  noise,  Figure  3.2  starts  by  comparing  the  CPM 
signal  and  its  noisy  version.  Figure  3.3  shows  the  performance  of  the  A4  estimator 
when  processing  this  noisy  signal.  What  makes  this  plot  suspect  is  that  increas¬ 
ing  the  number  of  symbols  does  not  register  a  cleaner  maximum.  To  illustrate  the 
performance  over  all  symbol  timings,  Figure  3.4  displays  the  function: 

1— >  Norm[A4(x|e)], 

where  the  normalization  scales  A4  to  the  interval  [0,1]: 

A  (  n  A4(x|e)  -  min{A4(x|g)} 

4X6  max{A4(x|e)}  —  min{A4(x|e)} 

An  ideal  estimator  would  peak  along  the  45°  reference  line  drawn  on  the  image. 
Figure  3.4  shows  that  A4  is  a  reasonable  estimator  for  this  large  SNR. 
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Scaled  A4(x|e) 


Maximum  Likelihood  Estimate  for  Symbol  Timing  t-0 


e:  /j=[4/16  5/16];  A$ps=10;  SNR=0  dB 


Figure  3.3:  Robustness  of  the  A4(x|e)  estimator. 
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3.1.2  A  Superbaud  Timing  Estimator 

The  maximum  likelihood  function  for  the  superbaud  timing  77= 0,  1,  . . . ,  H  —  1  is  [10, 
Eq.  4.1]: 

A„(x,  e\fj)  =  xH[F(k ,  -  k2)Ts,  k\Ts  -  (fj  +  e)T)]x, 

where  e  denotes  the  maximum  likelihood  estimate  for  e. 

Figure  3.5  displays  the  baseband  received  and  transmitted  CPM  signals  with  the 
timing  errors: 

r  =  (fj  +  e)T  fj  =  0;  e—  —0.3. 

Figure  3.5  shows  the  symbol  timing  estimator  e  is  tracking  e.  Using  this  estimate  e, 
Figure  3.5  presents  the  superbaud  timing  estimator  where  the  lines  are  drawn  between 
the  integer  values  of  ?/  for  clarity. 


CPM  Simulations 
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Figure  3.5:  CMP  signal  in  noise. 


24 


Maximum  Likelihood  Estimate  for  Symbol  Timing  f=-0.3;  eML=-0.3 


e:  /,=[ 4/16  5/161;  ^=10:  A/sym-5;  SNR=20  dB 

Figure  3.6:  Symbol  timing  estimate. 


Maximum  Likelihood  Estimate  for  Super  Baud  r|=0;  rjMI  =0 


n:  *=[4/16  5/16];  A^sps=10;  /Vsym=5;  SNR=20  dB;  e=-0.3;  eMl  =-0.3 


Figure  3.7:  Superbaud  estimate. 
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Encouraged  by  this  correct  result,  we  experimented  by  changing  the  superbaud 
timing  to  fj  =  1.  Figures  3.8,  3.9,  and  3.10  show  the  time  series,  the  failed  symbol 
time  estimate,  and  the  resulting  failure  of  the  superbaud  timing  estimate. 


CPM  Simulations 


Figure  3.8:  CMP  signal  in  noise;  fj  =  1. 
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448 


Maximum  Likelihood  Estimate  for  Symbol  Timing  e=-0.3;  eMI  =0.45 


Figure  3.9:  Symbol  timing  estimate;  fj  =  1. 
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It,  is  worth  asking — do  more  symbols  improve  this  estimator?  Instead  of  using 
only  5  symbols,  the  number  of  symbols  was  increased  to  20.  Figures  3.11,  3.12,  and 
3.13  show  the  time  series,  the  failed  symbol  time  estimate,  and  the  resulting  failure 
of  the  superbaud  timing  estimate.  This  result  may  explain  the  running  update  of  fj 
found  in  [10]. 


CPM  Simulations 


samples:  SNR=20  dB;  t=7  samples;  A?sps— 10 


Figure  3.11:  CMP  signal  in  noise;  ij  =  1;  20  symbols. 
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Maximum  Likelihood  Estimate  for  Symbol  Timing  e=-0.3;  eM[  =0.5 
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Figure  3.12:  Symbol  timing  estimate;  fj  =  1;  20  symbols. 


Maximum  Likelihood  Estimate  for  Super  Baud  r|=l ;  rjM,  =0 
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Figure  3.13:  Superbaud  estimate;  77  =  1;  20  symbols. 
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3.2  2 -h  Timing 

An  approach  specialized  t.o  2 -h  delay  is  the  maximum-likelihood  estimator  [29,  Eq. 
16]: 

N, 

AP(x|f)  :=  |AT|2;  X  :=  ^  xB(kTs)F(h0,  h\ ,  kTs  -  r), 

k= 1 

where  the  ubiquitous  F  function  is  [29,  Eq.  15]: 

Nsym/ 2—1 

F(ho,  h] ,  t)  :=  n  cos(277h0q(t  -  2 k.T))  cos(47r h0q(t.  -  2k:T)) 

k=0 

x  cos(2nhiq(t  —  (2k  +  1  )T))  cos(47rhi<3'(t  —  (2k  +  1)7")). 

The  product  indices  are  found  on  [29,  page  17].  The  “region  of  uncertainty”  for  t  is 
the  interval  [0,  T]  so  the  maximum-likelihood  estimator  [29,  page  5] : 

f  :=  argmax{Ap(x|f)  :f  G  [0,  T]} 

is  actually  an  estimate  of  the  symbol  timing  parameter  e  shifted  by  T/2. 
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The  robustness  of  f  is  problematic.  Figure  3.14  displays  a  CPM  signal  and  its 
noise-corrupted  version.  To  test  the  estimator  to  its  best  advantage,  the  delay  is 
restricted  to  r  €  [0,  T ]  and,  as  the  figure  shows,  the  noise  is  very  small. 


CPM  Simulations 


time:  SNR-20  dB;  t=0. 5;  N  =10;  TV  =100 

sps  ’  sym 

Figure  3.14:  Delayed  CPM  and  its  noisy  version  observed  over  100  symbols. 
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Figure  3.15  plots  AP(x|r)  over  a  symbol  interval  [0, 1].  The  actual  value  of  r  is 
marked  on  the  plot.  The  monotonic  shape  of  AP(x|r)  is  explained  by  opening  up  the 
plot. 
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Figure  3.15:  AP(x|r)  over  a  symbol  interval. 


Figure  3.16  plots  Ap(x|r)  over  the  hill  observation  interval.  The  lack  of  periodicity 
precludes  various  foldings  to  [0,  T]  to  rescue  this  estimator. 


Maximum  Likelihood  Estimate  for  x=0.5 


t:  /»=[ 4/16  5/16];  A^-10;  A^^lOO;  SNR-20  dB 

Figure  3.16:  Ap(x|r)  over  the  observation  interval. 
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3.3  Frequency  Synchronization 

The  maximum  likelihood  estimator  for  frequency  is  the  quadratic  form  [28,  Eq.  16]: 

Apw(x|/)  :=  (z  ©  x)T  H  (z  0  x). 

Here  x  is  the  vector  of  time  samples  of  the  noisy  CPM  signal 


x:=rB(t)  (t  :=  [kTs]). 


Recall  that  the  number  of  time  samples  TV,  is  the  product  of  the  number  of  observed 
symbols  and  the  number  of  samples  per  symbol:  TV,  =  TVsyrn  x  TVsps.  The  complex 
sinusoid  is 

z  :=  exp(-j27r/t). 

The  “©”  denotes  the  Hadamard  or  pointwise  vector  product.  The  matrix  H  is  TV,  x  TV, 
with  [28,  Eq.  17]: 

,  .  .  _  J_  f2T  ^  sin(8?r hop(k2Ts  -  t  -  2??T,  (k2  -  ki)Ts)) 

15  1  ’  2T  Jo  n^_ooAsm(27Th0p(k2Ts  ~  t  -  2nT,(k2  -  ki)Ts)) 

sin(87r liipjhTs  -t  —  2nT,  {kg  ~  ^i)^))  ^ 

4sin(27rhip(A:2T’s  —  t  —  2nT,  (&2  —  A:i)Ts))  ’ 


where  [28,  Eq.  18]: 

p{t,At.)  :=  q(t)  -  q{t  -  At). 

With  some  work,  one  can  verify  that  H  is  real,  symmetric,  Toeplitz,  and  positive. 
Observe  that  the  conjugation  is  twisted  in  this  form.  Taking  the  complex  conjugate 
converts  the  twisted  form  into  a  proper  quadratic  form.  However,  we  use  the  original 
form  in  the  simulations. 

The  following  figures  illustrate  this  maximum  likelihood  estimator.  The  received 
signal  is 

rB(t)  =  ej2*fo,sB(t  -  t)  +  g(t). 

As  in  the  preceding  simulations,  the  symbol  period  T  and  sampling  period  Ts  are  set 
as 

T  =  1,  Ts  =  0.1 

so  that  there  are  10  samples  per  symbol. 
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Figure  3.17  compares  the  transmitted  and  received  multi-/?,  signals.  The  upper 
and  lower  panels  display  the  real  (T)  and  imaginary  (Tv)  parts  of  these  baseband 
signals.  The  noise  is  significant,  and  the  number  of  symbols  is  relatively  small. 


CPM  Simulations 


time:  SNR=5  dB;/o=0.2;  x=0.38;  Wsym=20 

Figure  3.17:  Comparing  the  multi-/?  signal  sB(t)  and  its  noisy,  frequency-shifted 
version  rB(t). 
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Although  Figure  3.17  shows  the  received  signal  is  corrupted  by  noise  and 

shifted  in  frequency,  the  maximum  likelihood  estimate  shown  in  Figure  3.18  clearly 
registers  a  single  maximum  that  is  relatively  close  to  the  frequency  offset.  This 
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Figure  3.18:  A  maximum  likelihood  estimate  for  the  frequency  offset. 

plot  indicates  that  this  frequency  estimator  is  robust,  accurate,  and  requires  only  a 
relatively  small  number  of  symbols. 


Maximum  Likelihood  Estimate  for  7^=0. 24*:  ./mle~0-  15 1 52T 

- T - , - , - - T - 

T 


I _ i _ i _ i _ i _ i _ i - 1 — — - j — ■ — — — i — ■  — — J 

1  -0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1 

/  h=[ 4/16  5/16];  t=0.38;  SNR=5  dB;  N  =10;  JVsym=20 
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To  explore  this  claim,  Figure  3.19  plots  the  maximum  likelihood  estimate  /ml  as 
a  function  of  the  frequency  offset.  The  solid  line  on  this  plot  displays  the  target,  of  the 


Maximum  Likelihood  Estimate  for  frequency  offset 


»  1 - - -  *■ - 1 - L - 1 - -L— - 1  - -  1 

-0.5  -0.4  -0.3  -0.2  -0.1  0  0.1  0.2  0.3  0.4  0.5 

f-h=[ 4/16  5/16];  t=0.38;  SNR=20  dB;  N= 20 ;  /V  =10 

*'0  L  9  sym  ’  sps 


Figure  3.19:  A  maximum  likelihood  estimate  of  the  frequency  offset — high  SNR. 

/mle’s.  Over  each  frequency  offset,  the  column  of  dots  are  the  /mle’s  generated  for 
several  trials.  The  “grid”  pattern  is  visible  because  the  /mle’s  are  generated  on  the 
same  grid  as  the  frequency  offsets.  The  variance  in  the  /mle’s  appears  to  be  relatively 
constant  and  small  over  the  frequency  range.  Thus,  the  claims  of  an  accurate  and 
robust  estimator  are  validated — at  this  high  SNR. 
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Figure  3.20  runs  the  same  simulations  using  a  low  SNR.  As  before,  the  variance 
in  the  /mle’s  appears  to  be  relatively  constant  over  the  frequency  range.  Thus, 
this  maximum  likelihood  estimate  of  the  frequency  offset  appears  to  be  robust  and 
accurate  on  relatively  short  data  segments. 


Maximum  Likelihood  Estimate  for  frequency  offset 


fo:h=[ 4/16  5/16];  x=0.38;  SNR-5  dB;  Nsym=2Q;  =10 


Figure  3.20:  A  maximum  likelihood  estimate  of  the  frequency  offset — low  SNR. 
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3.4  Discussion 


CPM  synchronizers  have  been  under  development  for  decades  and  a  considerable 
literature  is  available.  A  standard  approach  is  to  start  from  the  joint  maximum 
likelihood  estimator  and  factor  the  joint  estimator  into  several  single-parameter  es¬ 
timators  using  independence  assumptions.  The  preceding  synchronizers  all  estimate 
a  single  parameter.  For  future  efforts,  Table  3.1  highlights  a  few  of  the  joint  estima¬ 
tors.  Synchronization  can  also  key  off  a  preamble  or  specified  data  sequence.  Such 


Table  3.1:  CPM  Synchronization  Literature.  NDA=Nondata  Aided.  ML=maximum 
likelihood. 


Multi-# 

NDA 

Channel 

Parameters 

Algorithm 

Reference 

V 

AWGN 

delay-)- phase 

ML 

[6],  [23] 

V 

v/ 

AWGN 

phase+data 

[19] 

x/ 

x/ 

AWGN 

delay+phase+data 

[31] 

v/ 

x/ 

AWGN 

delay+phase 

ML 

[15] 

x/ 

AWGN 

frequency,  delay 

[42] 

x/ 

AWGN 

delay+phase 

ML/Walsh 

[39] 

x/ 

AWGN 

frequency 

[4] 

synchronizers  are  called  data-aided.  The  synchronizers  covered  in  this  report  are 
nondata-aided.  As  the  table  shows,  the  synchronization  literature  is  focused  on  the 
additive  noise  channels. 
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4 

Synchronizer  Performance 


This  section  generalizes  the  simple  CPM  receiver  (CMPO)  of  Section  2  to  handle 
timing,  phase,  and  frequency.  Section  4.1  extends  CPMO  to  handle  the  phase  and 
timing  (CPM1).  Section  4.2  reports  on  the  performance  of  CPM1,  and  we  conclude 
that  phase  and  timing  for  multi-/?  CPM  in  Gaussian  noise  are  readily  handled  with 
this  simple  design  of  CPM1.  Thus,  a  more  complex  CPM  receiver  must  beat  CPM1 
by  several  dB  to  justify  trading  complexity  for  performance. 

Section  4.3  considers  several  approaches  for  frequency  synchronization.  One  sim¬ 
ple  approach  of  estimating  the  frequency  offset  is  developed  in  Section  4.4.  This 
generalizes  CPM1  to  handle  timing,  phase,  and  frequency  shifts  (CPM2).  Section  4.5 
illustrates  the  performance  and  limitations  of  CPM2.  If  too  few  symbols  are  ob¬ 
served,  an  inaccurate  estimate  of  the  frequency  shift  forces  a  noise  floor  in  CPM2. 
This  observation  raises  excellent  engineering  questions: 

•  How  sensitive  is  a  CPM  receiver  to  frequency  errors?  That  is,  how  good  does 
the  frequency  synchronizer  have  to  be? 

•  What  is  the  tradeoff  between  frequency  synchronization  and  the  error  rate? 
That  is,  given  a  desired  error  rate,  how  accurate  does  the  frequency  synchronizer 
have  to  be? 

Section  4.6  closes  with  a  discussion  of  these  questions  and  the  usage  of  these  CPM 
receivers  to  assess  synchronizer  performance. 


4.1  CPM1:  A  Simple  CPM  Receiver 

Section  2  compared  the  performance  of  the  simple  CPM  receiver  (CPMO)  against 
designs  in  the  IEEE  literature — assuming  perfect  synchronization.  The  performance 
of  CPMO  was  credible.  This  section  generalizes  CPMO  to  CPM1  that  handles  timing 
r  and  phase  errors  Ac/).  Credible  performance  is  reported  in  the  following  section. 
Frequency  errors  are  treated  later  so  we  assume  A/c  =  0  for  CPM1. 
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Assume  that.  Nsym  symbols  have  been  transmitted: 

^sym 

sB(t)  =  exp Mt)  =  2tt  h[k]akq{t  -  kT ). 

k=l 

Assume  the  received  baseband  signal  rB(t)  is  corrupted  by  additive  noise  and  demod¬ 
ulation  errors  of  the  following  form: 

r e(t)  =  -  r)  +  gB(t). 

The  received  phase  then  has  the  following  form: 

^R{f)  =  A0  +  27tA fct  +  -  r)  + 

where  the  last  term  is  the  phase  noise  induced  by  the  additive  noise  {<?e(t)}.  The 
problem  is  to  recover  and  synchronize  i(^B(t)  from  rB(t). 

To  introduce  the  CPM1  design,  recall  that  £{t)  denotes  the  piecewise  linear  func¬ 
tion 

OO 

m  ■=  2tt  £  ~  kT)2T- 

k~ — oo 

With  L  =  1,  the  first,  derivative  is 

OO 

f(t)  ■■=  2*  E  A  W<  -  kT). 

k~— oo 

Recall  also  that  C  denotes  collection  of  all  these  piecewise  linear  functions.  The  CPM1 
design  is  a  generalization  of  the  CPMO  design  where  £(t  -  r)  is  fit  to  the  phase  i>R(t) 
so  the  instantaneous  frequencies  fk  s  can  be  read  off  from 

CPM1-1  Estimate  the  phase  of  sB(t) 

1.  Phase  unwrap  [20]: 

:=  unwrap(angle(rs(t))). 

2.  Compute  the  instantaneous  frequency  to  remove  the  phase  error: 

i,'R{t)  —  2ttA  fc  +  ij/B(t  —  t)  + 

By  assumption,  A/c  =  0.  If  the  additive  noise  is  small, 

V4W  &ip'B(t-T). 

3.  Integrate  to  get  an  estimate  of  ijjB{t  —  r): 

i>B(t  -  t)  :=  f  . 

Jo 
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CPM1-2  Synchronize: 

1.  Estimate  the  symbol  timing  error  as 

{e,  1}  :=  argmin{||^B(°  ~  T)  ~  £(°  ~  :  f  £  [  1/2, 1/2],  f  €  £}• 

That  is,  f  is  a  symbol  time  that  best  explains  the  phase  as  a  piecewise 
linear  function. 

2.  Estimate  the  instantaneous  frequency  from  the  best  piecewise  linear  min- 
imizer  (L  =  1): 

fB(t)  :=  2-?(«  -  ST). 

3.  Estimate  the  superbaud  error  rj:  If  the  symbol  timing  is  correct, 

.fk  ~ 

so  choose  the  delay  that  minimizes 


'  3  ' 

(  Nsym  | 

1 

fj  :=  argmin  <  II  fk/h+v  -  a||_oo  :  r]  =  0, 1  >  , 

a  = 

-1 

[  k=  1  J 

.  -3  . 

CPMl-3  Decode: 

ak  :=  argmin{|  fk/hk+~  -  a\  :  a  =  ±1,  ±3} 


42 


For  orientation,  the  following  figures  illustrate  CPM1  in  a  typical  simulation. 
Figure  4.1  compares  the  transmitted  (no  noise)  CPM  signal  sB(t)  against  the  noisy, 
time-shifted,  and  phase-rotated  baseband  received  signal  Although  the  noise 

is  still  relatively  small,  the  phase  does  distort  the  received  signal. 
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Figure  4.1:  Transmitted  and  received  CPM1  signals. 
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Figure  4.2  makes  the  corresponding  comparison  on  the  phases  and  instantaneous 
frequencies.  Even  though  the  SNR  is  large,  and  phases  are  close,  the  derivative 
amplifies  the  phase  noise  in  the  instantaneous  frequency. 


CPM  Timing  Only:  phase  &  frequency  estimates 


X=3.8  samples;  A<|>=-60  deg 


samples:  SNR=10  dB;  Nsp  =10;  Vsym=20 


Figure  4.2:  Transmitted  and  estimated  CPM1  instantaneous  phase  and  frequency. 
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Figure  4.3  shows  that  the  symbol  timing,  the  superbaud  timing,  and  the  decod¬ 
ing  has  been  successful  despite  the  noise  in  the  instantaneous  frequency.  The  step 
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Figure  4.3:  CPM1  performance. 

function  (blue  for  color  viewers)  plots  the  instantaneous  frequency  of  the  transmit¬ 
ted  signal:  The  jagged  (green)  line  is  the  received  phase  after  correcting 

for  phase,  frequency,  and  superbaud  timing.  The  (blue)  diamonds  mark  the  instanta¬ 
neous  frequency  estimated  from  the  received  phase.  Figure  4.3  reports  that  estimated 
probability  of  symbol  error  is  zero:  Pe  =  0.  From  that  observation,  we  infer  that  the 
estimate  of  the  symbol  timing  did  also  work. 
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4.2  CPM1:  Timing  and  Phase  Performance 

The  performance  of  CPM1  is  illustrated  in  the  following  plots.  Figures  4.4  and  4.5 
plot  estimates  of  the  probability  of  symbol  error  as  a  function  of  the  SNR  and  Ei,/Nq, 
respectively.  Equation  2.2  is  used  to  map  the  SNR  to  Eb/N0.  In  both  plots,  the  timing 


CPM  Timing  Only:  Probability  of  Symbol  Error 


Figure  4.4:  CPM1  SNR  Performance;  timing  errors  only. 

error  is  held  fixed  at  r  =  3.8  samples  and  the  CPM  simulations  contain  Nsym  =  100 
symbols.  Both  plots  present  each  simulation  as  a  dot  and  draw  the  mean  probability 
of  error  as  line.  These  plots  demonstrate  that  CPM1  is  at  least  functional  with  respect 
to  timing.  A  nice  feature  of  the  CPM1  design  is  that  the  phase  error  A <f>  is  eliminated 
by  the  derivative.  Figures  4.6  and  4.7  verify  that  CPM1  functions  regardless  of  the 
phase  error. 
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Peh=[4!  16  5/16] 


CPM  Timing  &  Phase  Only:  Probability  of  Symbol  Error;  t=3.8  samples;  A ({>—60  deg 


SNR  (dB):  ^,,=  1000;  /Vsym=100 


Figure  4.6:  CP  Ml  SNR  Performance;  timing  and  phase  errors  only. 


CPM  Timing  &  Phase  Only:  Probability  of  Symbol  Error;  r=3.8  samples;  A<{>=-60  deg 


V"o  <dBI;  ^„p-i°00;  V100 


Figure  4.7:  CPM1  E^/Nq  Performance;  timing  and  phase  errors  only. 
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4.3  Designs  for  Timing,  Phase,  and  Frequency 

The  reason  that,  CPM1  handled  the  phase  error  A<p  is  that  taking  the  derivative  of 
the  received  phase 

V4C0  =  V4(*  -  T)  +  g(* ) 

automatically  removes  the  phase  error — assuming  that  no  frequency  error  A fc  is 
present.  Recall  that  tpG^)  denotes  the  phase  error  induced  by  the  additive  noise. 
If  the  frequency  error  is  present,  the  instantaneous  frequency  of  the  received  CPM 
signal  rs{t)  is 

V4W  =  2ttA/c  +  xJ/B(t  -  r)  + 

Thus,  it  is  natural  to  take  another  derivative  to  remove  the  frequency  error: 

V4  (<)  =  V4(*  -  T)  +  V4(*)- 

The  problem  is  that  using  the  1REC  phase  pulse, 

Jit)  =  ^rMo,T\(t)  ==>  q'\t)  =  -  6{t  -  T)}  . 

That  is,  most  of  the  signal  is  also  wiped  out  and  the  information  is  carried  in  two 
noise-corrupted  estimates  of  the  delta  functions  that  now  overlap  with  the  preceding 
and  following  pulse. 

Simulations  verified  that  using  two  derivatives  with  the  1REC  pulse  did  compen¬ 
sate  for  the  frequency  offset — but  the  symbol  error  rate  was  too  large.  Nevertheless, 
several  schemes  to  compensate  for  the  frequency  shift  come  to  mind. 

•  Instead  of  using  the  1REC  phase,  use  a  phase  pulse  with  a  non-vanishing  second 
derivative  (i.e.,  the  chirps  in  [7,  Fig.  2.4]). 

•  Estimate  A  fc  from  the  sample  mean  of 

•  Simultaneously  estimate  A fc  and  the  /it’s  using  a  Viterbi  decoder. 

•  Use  the  frequency  estimator  of  Section  3.3. 

Because  we  are  looking  for  the  simplest  design,  the  next  section  explores  the  CPM 
receiver  that  uses  the  sample  mean. 
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4.4  CPM2:  Timing,  Phase,  and  Frequency 

CPM2  generalizes  CPM1  to  handle  frequency  offsets  using  the  sample  mean  of  the 
received  phase.  The  received  signal 

rB{t)  =  ejA<t,ej2nAfc1sB(t  -  r)  +  gB(t) 

has  instantaneous  frequency 

=  2vrA/r  +  V4(*  -  r)  +  V4(0 


with  expected  value 

£%(*)]  =  2ttA/c  +  £%(*  -  r)]  +  EWc{t)]  =  2ttA/c, 

using  the  circularity  of  the  noise  and  assuming  all  symbols  a*,  =  ±1,  ±3  are  equally 
likely.  If  a  suffi.cie.ntly  large  number  of  symbols  are  observed,  a  sufficiently  accurate 
estimate  of  the  frequency  offset  is  obtained  from  the  sample  mean.  The  CPM2  design 
is  a  generalization  of  the  CPM1  design  using  this  estimate  of  the  frequency  offset. 

CPM2-1  Estimate  the  phase  of  sB(t.) 

1.  Phase  unwrap  [20]: 


Mt)  :=  unwrap(angl e(rB(t))). 

2.  Compute  the  instantaneous  frequency  (IF)  to  remove  the  phase  error: 

xj/R{t)  =  2ttA/c  +  -t)  +  il>G(t). 

3.  Estimate  the  frequency  error: 

A/c  :=  ^-mean [rf/R(t)]. 


4.  Integrate  to  get  a  clean  phase: 

Mt  ~  t)  :=  ^  {V4(0  “  27tA/c}  dt'. 


CPM2-2  Synchronize: 

1.  Estimate  the  symbol  timing  error  as 

{ej}  :=  argmin{||^B(o  -  r)  -  £(o  -  eT) ||2  :  e  €  [-1/2, 1/2],  ^  6  £}. 
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2.  Estimate  the  instantaneous  frequency  from  the  best  piecewise  linear  min- 
imizer  (L  =  1): 

M<)  “  2-?(«  -  KT). 

3.  Estimate  the  superbaud  error  rj:  If  the  symbol  timing  is  correct, 


fk  ~  O'^.  /?,[£. ^_r^| 

so  choose  the  delay  that  minimizes 


f]  :=  argmin 


!^sym 

Z  II  fk/hk+r,  *^||  —  oo 
k=l 


:  V  =  0, 1 


CPM2-3  Decode: 

ak  :=  argmin{|/fc//?*+~  -  o  |  :  a  =  ±1,  ±3} 
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The  following  figures  illustrate  CPM2  in  a  typical  simulation.  Figure  4.8  compares 
the  transmitted  (no  noise)  CPM  signal  sB(t)  against  the  noisy  time-shifted,  phase- 
shifted,  and  frequency-shifted  baseband  received  signal  rB(t).  The  noise  is  small  but 
the  frequency  offset  is  spinning  the  received  signal. 

CPM2  Baseband  Simulations 
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Figure  4.8:  Transmitted  and  received  CPM  signals. 
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Figure  4.9  makes  the  corresponding  comparison  on  the  phases  and  instantaneous 
frequencies.  Even  though  the  SNR  is  large,  and  phases  are  close,  the  derivative 
amplifies  the  phase  noise  in  the  instantaneous  frequency. 


CPM2:  phase  and  frequency  estimates 
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Figure  4.9:  Transmitted  and  CPM2  estimate  of  the  instantaneous  phase  and  instan¬ 
taneous  frequency. 
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Figure  4.10  compares  the  spectra  of  the  transmitted  CPM  signal  the  re¬ 

ceived  signal  i'[j  (t).  and  the  estimated  signal 

sE(t)  ■=  exp -  r)). 

The  spectrum  of  received  signal  is  registering  the  frequency  shift.  The  frequency 
estimate  appears  to  he  close  to  the  frequency  shift.  We  see  that  the  spectrum  of  the 
estimated  signal  overlaps  the  spectrum  of  the  transmitted  signal. 


frequency  (Hz):  SNR=20  dB;  Vsp  =10;  Vsym=20;  A/Est=0.97695  Hz 

Figure  4.10:  CPM2  spectra. 
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Figure  4.11  shows. that  the  symbol  timing,  the  superbaud  timing,  and  the  decod¬ 
ing  has  been  successful  despite  the  noise  in  the  instantaneous  frequency.  The  step 


Instantanous  frequencies:  A/^l  Hz;  t=3.8  samples;  A(|>=-60  deg 
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Figure  4.11:  CPM2  performance. 

function  (blue  line  segments)  shows  the  instantaneous  frequency  of  the  transmitted 
signal:  The  noisy  (green)  line  is  the  received  phase  after  correcting  for  phase, 

frequency,  and  superbaud  timing.  The  (blue)  diamonds  mark  the  instantaneous  fre¬ 
quency  estimate  of  that  estimated  symbol  period  obtained  from  the  noisy  received 
phase.  That  the  estimated  probability  of  symbol  timing  is  zero  (Pe  =  0)  means  that 
the  symbol  timing  functioned. 
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4.5  CPM2:  Timing,  Phase,  and  Frequency  Performance 

Figure  4.12  reports  on  the  performance  of  CPM2.  As  in  the  preceding  simulations, 
the  delay  is  3.8  sample  periods  (r  =  3.8  x  Ts )  and  the  phase  shift  is  A (j)  =  —60°. 
The  frequency  shift  is  scaled  to  the  sample  period  of  Ts  =  0.1.  Using  seconds,  we 
shift  the  CPM  signal  by  A /  =  1  Hertz.  The  effect  of  the  frequency  shift  on  the 


()  Probability  of  Symbol  Error  for  CPM2:  A f=  1  Hz:  r=3.8  samples:  A0=-6O  deg 
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SNR  (dB):  A^rep-2000:  ^=100 


Figure  4.12:  CPM2  performance  using  only  100  symbols. 

simple  synchronizer  is  clear:  As  the  noise  decreases,  a  noise  floor  is  revealed.  Thus, 
the  frequency  estimate  is  too  coarse  for  a  credible  receiver.  However,  the  frequency 
estimator  was  gotten  from  100  symbols.  Would  taking  more  symbols  improve  the 
performance  of  CMP2?  Figure  4.13  shows  collecting  200  symbols  does  improve  the 
performance  of  CPM2.  Doubling  the  number  of  symbols  pushed  the  error  floor  down 
10  dB.  However,  an  error  floor  still  exists  and  is  still  too  high.  Figure  4.14  reports  that 
opening  the  observation  interval  to  500  symbols  continues  to  improve  the  performance 
of  CPM2. 

Before  we  leave  this  section,  it  is  worthwhile  to  put  these  frequency  shifts  in 
context.  A  frequency  shift  of  A/  =  1  Hertz  is  commensurate  with  the  symbol  period 
T  =  1  second.  Scaled  to  a  real-world  system  using  T  =  1  //s,  the  frequency  shift  is 
A/  =  1  MHz.  This  exceeds  typical  demodulation  and  frequency  errors  by  several 
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Probability  of  Symbol  Error  for  CPM2:  A/M  Hz:  x=3.8  samples:  A<|>=-60  deg 


SNR  (dB):  10000;  7Vsyni=200 

Figure  4.13:  CPM2  performance  using  200  symbols, 
orders  of  magnitude  [2] . 

4.6  Discussion 

The  progression  of  CPM0,  CMP1,  and  CPM2  shows  how  to  generalize  the  very  simple 
multi-/?,  CPM  receiver  to  handle  timing,  phase,  and  frequency  synchronization.  The 
progression  also  shows  that  the  synchronizer’s  performance  should  be  measured  in 
the  context  of  a  receiver.  For  example,  increasing  the  accuracy  of  the  frequency  syn¬ 
chronizer  by  10%  is  meaningless  until  the  10%  improvement  maps  to  an  improvement 
in  its  CPM  receiver. 

The  simplicity  of  these  CPM  designs  also  lets  us  swap  out  various  synchronizers. 
For  example,  the  frequency  synchronizer  of  Section  3.3  can  estimate  the  frequency 
A/e  used  in  CPM2-1.  We  can  ask:  Does  the  accuracy  of  the  more  complex  estimator 
let  us  shorten  the  observation  interval?  More  generally,  we  could  omit  a  frequency 
estimator  entirely.  Indeed,  we  could  simply  hand  the  CPM  receiver  the  true  frequency 
A  fc  corrupted  by  Gaussian  noise: 

A/c  =  A/c  +  A/. 
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Probability  of  Symbol  Error  for  CPM2:  zW  I  Hz:  x=3.8  samples;  A<|) =-f>0  deg 


SNR  (dB):  /Vrep=  10000;  Wsym=500 

Figure  4.14:  CPM2  performance  using  500  symbols. 

The  performance  of  the  CPM  receiver  could  then  be  known  as  function  of  the  vari¬ 
ance  of  A/.  Knowing  this  function  sets  the  accuracy  requirements  of  any  frequency 
estimator.  However,  tinkering  with  frequency  estimators  in  additive  noise  is  much 
less  of  a  problem  than  multipath  channels. 
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5 

Multi-h  in  Multipath 


The  received  signal  in  multipath  is  a  sum  of  delayed  and  faded  versions  of  the  trans¬ 
mitted  signal: 

M 

rB{t)  =  a™.(t)sB(t  -  Tm) 

777=1 

Each  fading  process  {om(t)}  modulates  the  CPM  signal.  The  fading  processes  {«„,.(?)} 
under  consideration  have  power  levels  ranging  from  0  to  -10  dB.  We  have  seen  that 
the  CPM  receivers  need  a  signal  in  excess  of  10  dB  over  additive  noise.  Consequently, 
some  fading  processes  may  inject  too  much  noise  for  these  simple  CPM  receivers. 
Sections  5.1  and  5.2  verify  this  noise  limitation  using  a  UHF  channel  model  of  the  San 
Diego  Harbor.  Good  exploitation  of  multipath  may  require  simultaneously  estimating 
both  the  channel  and  signal  [2].  Section  5.3  formalizes  this  observation  and  concludes 
this  report. 

5.1  Harbor  Multipath 

A  communications  channel  can  be  modeled  as  a  stochastic,  time-varying  linear  filter 
h(t,  t)  that  maps  the  baseband  transmitted  signal  sB(t)  to  the  received  signal  rB(t) 
as  [9]: 

/OO 

h(t,t-  t)sb(t)cIt. 

-OO 

A  standard  multipath  model  is  the  Quadrature  Modulation  Fading  Simulator  (QMFS) 

[32]: 

M 

h(t,r )  =  ^  a,m{t)6{r  -  rm). 

77?  =  1 

The  received  signal  is  a  sum  of  delayed  and  faded  versions  of  the  transmitted  signal: 

M 

^*b(0  ^"777.) 

777,= 1 
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Typically,  the  fading  processes  {«,„(/)}  are  complex- valued,  jointly  wide-sense  sta¬ 
tionary  (JWSS),  narrow-band,  Gaussian  random  processes: 

am(t)  =  ePK'Dt{'rm  +  gm(t)}. 

The  Dopper  shift  is  denoted  by  //>.  The  specular  component  models  coherent  re¬ 
flection  using  the  complex- valued  reflection  coefficient  7m.  The  diffuse  component 
is  modeled  by  the  narrow-band,  zero- mean,  Gaussian  random  process  {#„,(/)}  with 
variance  <r2.  The  K  factor  for  the  rath  path  is 


Random  phases  are  simultaneously  absorbed  in  jm  and  {<7m  (/•)}•  For  the  delays, 
general  models  take  the  rm’s  as  time- varying.  Constant  delays  sufficed  for  modeling 
the  channels  in  the  Extended  Littoral  Battlefield  (ELB)  [2].  The  delays  and  statistics 
of  the  fading  processes  determine  the  channel  model.  Table  5.1  lists  the  statistics  for 
the  Two-Path  Harbor  Model.  All  the  measurements  are  referenced  to  Path  1. 

Table  5.1:  Two-Path  Harbor  Model  [3]. 


Path 

Delay 

(IS 

Power 

dB 

Doppler  shift 

Hz 

Fade  rate 

Hz 

Fading 

1 

0.0 

0 

0 

1.6 

Rican  K  =  136 

2 

2.0 

-5 

-5 

3.8 

Rican  K  —  14 

5.2  Multi-/i  in  the  Harbor  Multipath 

This  section  pushes  the  multi-//  CPM  signal  through  the  Harbor  multipath  of  Ta¬ 
ble  5.1  and  then  attempts  to  processes  the  received  signal.  To  match  the  bandwidth 
of  the  ELB  measurements,  106  symbols  per  second  are  transmitted.  The  symbol 
interval  is  then 

T  =  1  (/is) 

If  the  observation  interval  is  100  symbols,  the  large  K  factors  and  the  slow  fade  rates 
make  the  Harbor  multipath  almost  time  invariant: 

rB(t)  =  ai(t)sB(t  -  ti)  +  tt2(t)sB(t  -  t2) 

=  {7i  +01MWO  +  e-j2KfDt{j2  +  92(t)}sB(t  -  t2) 

~  7isr?(/)  +  e_j27r/D,72Sf?(/-r2). 
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Because  the  observation  interval  is  only  100//S,  the  5-Hertz  Doppler  is  absorbed  into 
72-  The  Two-Path  Harbor  Model,  assuming  100  symbols,  then  admits  the  approxi¬ 
mation 

rB(t)  =  7i sB(t)  +  72 sB{t  -  r2). 

With  a  2-symbol  delay,  Path  2  looks  like  5  dB  noise  to  Path  1.  The  preceding  section 
demonstrated  that  5  dB  noise  will  kill  the  CPM  receiver  by  forcing  a  probability  of 
symbol  error  Pe  of  approximately  15%.  Path  2  forces  a  similar  effect. 

Figure  5.1  shows  the  transmitted  and  multipath  received  signals.  Even  though 
the  additive  noise  is  20  dB  down  from  Path  1,  Path  2  is  only  5  dB  down.  Compared 
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Figure  5.1:  CPM  in  Harbor  multipath. 

to  the  low-noise  plots  of  Section  4,  Path  2  is  a  significant  noise  term.  Figure  5.2 
compares  instantaneous  frequency  of  the  transmitted  signal  to  the  multipath  version. 
Path  2  has  injected  large  spikes  that  degrade  the  CPM  receiver  (CMP2)  and  returns 
a  probability  of  symbol  error  commensurate  with  5  dB  noise.  Additional  simulations 
verify  that  the  power  in  Path  2  is  equivalent  to  adding  noise  at  that  power.  In  contrast 
to  the  multipath  that  actually  improved  the  performance  of  an  LPI  receiver  [3],  the 
2-symbol  delay  in  Path  2  degrades  the  CPM  receiver.  Consequently,  exploitation  of 
multipath  may  require  the  channel  and  signal  estimation  to  use  these  CPM  receivers 
[2]- 
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1.5 


CPM2  in  2-Path  Harbor  Model:  Instantanous  frequencies 
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Figure  5.2:  CPM2  performance  in  Harbor  multipath. 

5.3  Summary 

This  report  sets  out.  a  simulation  framework  that  allows  phase,  timing,  and  frequency 
synchronizers  to  be  swapped  in  and  out  of  a  CPM  receiver  to  assess  end-to-end  per¬ 
formance.  The  additive  noise  simulations  demonstrated  the  functionality  of  the  basic 
CPM  receivers  of  this  report.  However,  the  preceding  multipath  simulation  demon¬ 
strated  these  CPM  receivers  require  additional  processing  to  function  in  multipath. 
Table  5.2  organizes  selected  papers  on  CPM  performance  over  RF  channels.  The  table 
shows  few  papers  consider  multi-/?  CPM  in  the  full-up  multipath.  Consequently,  the 
simulation  efforts  of  Phase  3  should  focus  on  extending  the  CPM  receiver  to  handle 
multi-/?,  in  multipath.  The  simplest  approach  is  to  simultaneously  estimate  both  the 
CPM  signal  and  channel  [2].  The  received  and  demodulated  CPM  signal  can  then  be 
processed  by  these  basic  CPM  receivers  The  research  efforts  of  Phase  3  can  attack 
the  following: 

Q-l  How  many  symbols  are  needed  for  synchronization? 

Q-2  How  does  synchronization  depend  on  the  precision  of  the  estimators? 

Q-3  Does  there  exist  a  general  lower  bound  on  the  number  of  symbols  to  synchronize? 
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Table  5.2:  Multipath  CPM  Literature;  NDA:=Nondata-Aided;  DCD:=Data  and 
Channel  Demodulation;  TDL:=Tapped  Delay  Line;  DFSE:=Decision  Feedback  Se¬ 
quence  Equalization. 


Multi-/? 

NDA 

Channel 

Algorithm 

Reference 

V 

QMFS 

DCD 

[48] 

V 

TDL 

MLSE 

[47],  [13] 

V 

Rican 

Trellis 

[1] 

V 

TDL 

T 

[16] 

V 

TDL 

DFSE 

[13] 

y/ 

V 

TDL 

MLSD 

[17] 

V 

Rayleigh+CCI 

1-DPD 

[18] 

V 

TDL 

unwrapping 

[26] 

V 

V 

MILSATCOM 

DFSE 

[27] 

V 

V 

Satellite-Mobile 

MLSE 

[33] 

V 

Rayleigh 

Trellis 

[25] 

V 

V 

CCI 

MLSE 

[45] 

V 

V 

Rican-Mobile 

MLSE 

[46] 

Q-4  When  does  multipath  destroy  synchronization? 

Q-5  Would  demodulation  of  both  the  channel  and  signal  improve  CPM  performance 
in  multipath? 

Q-6  Can  a  hybrid  CPM  system  overcome  multipath? 

Q-l  suggests  that  the  CPM  receiver  should  be  made  adaptive  to  the  channel.  The 
question  implicitly  asks  how  the  performance  depends  on  the  SNR.  Q-2  brings  up 
the  issue  of  the  intertwining  to  the  various  synchronizers  in  the  full  CPM  radio.  For 
example,  it  makes  little  sense  to  improve  the  symbol  timing  estimator  if  the  CPM 
receiver  is  limited  by  super  baud  estimator.  Q-3  seeks  a  general  design  bound  to 
benchmark  any  synchronizer.  Q-4  is  answered  by  the  Harbor  multipath  example  in 
preceding  section.  The  final  questions  are  points-of-departure  for  Phase  3.  Q-5  is 
a  straight-forward  question  based  on  the  demodulation  results  of  [2],  Q-6  is  more 
ambitious  and  based  on  the  mix  of  CPM  and  the  excellent  multipath  properties  of 
OFDM  [40], 
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